Mon. Not. R. Astron. Soc. 000, () Printed 1 February 2008 (MN IAT^X style file vl.4) 



The kinematics and the origin of the ionized gas in 
NGC 4036 



P. Cinzano 1 , H.-W. Rix 2 ' 3 , M. Sarzi 1 , E.M. Corsini 4 , W.W. Zeilinger 5 
and F. Bertola 1 

1 Dipartimento di Astronomia, Universita di Padova, vicolo dell'Osservatorio 5, 1-35122 Padova, Italy 

2 Max Planck Institut fur Astrophysik, Karl Schwarzschild Strafie 1, D-85748 Garching bei Miinchen, Germany 

3 Steward Observatory, University of Arizona, Tucson, AZ-85721, USA 

4 Osservatorio Astrofisico di Asiago, Dipartimento di Astronomia, Universita di Padova, via dell'Osservatorio 8, 1-36012 Asiago, Italy 

5 Institut fur Astronomic, Universitdt Wien, Tiirkenschanzstrafie 17, A-1180 Wien, Austria 



Accepted .... Received in original form 



ABSTRACT 

We present the kinematics and photometry of the stars and of the ionized gas near 
the centre of the SO galaxy NGC4036. Dynamical models based on the Jeans Equation 
have been constructed from the stellar data to determine the gravitational potential 
in which the ionized gas is expected to orbit. Inside 10", the observed gas rotation 
curve falls well short of the predicted circular velocity. Over a comparable radial region 
the observed gas velocity dispersion is far higher than the one expected from thermal 
motions or small scale turbulence, corroborating that the gas cannot be following 
the streamlines of nearly closed orbits. We explore several avenues to understand 
the dynamical state of the gas: (1) We treat the g collisionless ensemble of 

cloudlets and apply the Jeans Equation to it; this modeling shows that inside 4" the 
gas velocity dispersion is just high enough to explain quantitatively the absence of 
rotation. (2) Alternatively, we explore whether the gas may arise from the 'just shed' 
mass-loss envelopes of the bulge stars, in which case their kinematics should simply 
mimic that of the stars, he latter approach matches the data better than (1), but still 
fails to explain the low velocity dispersion and slow rotation velocity of the gas for 
5" < r < 10". (3) Finally, we explore, whether drag forces on the ionized gas may aid 
in explaining its peculiar kinematics. While all these approaches provide a much better 
description of the data than cold gas on closed orbits, we do not yet have a definitive 
model to describe the observed gas kinematics at all radii. We outline observational 
tests to understand the enigmatic nature of the ionized gas. 

Key words: galaxies: elliptical and lenticular - galaxies: individual: NGC 4036 - 
galaxies: ISM - galaxies: kinematics and dynamics - galaxies: structure 



1 INTRODUCTION 

Stars and ionized gas provide independent probes of the 
mass distribution in a galaxy. The comparison between their 
kinematics is particularly important in dynamically hot sys- 
tems (i.e. whose projected velocity dispersion is comparable 
to rotation). In fact in elliptical galaxies and bulges the am- 
biguities about orbital anisotropies can lead to considerable 
uncertainties in the mass modeling (e.g. Binney & Mamon 
1982; Rix et al. 1997). 

The mass distributions inferred from stellar and gaseous 
kinematics are usually in good agreement for discs (where 
both tracers can be considered on nearly circular orbits), 
but often appear discrepant for bulges (e.g. Fillmore, Boro- 



son & Dressier 1986; Kent 1988; Kormendy & Westpfahl 
1989; Bertola et al. 1995b). There are several possibilities to 
explain these discrepant mass estimates in galactic bulges: 

(i) If bulges have a certain degree of triaxiality, depend- 
ing on the viewing angle the gas on closed orbits can either 
move faster or slower than in the 'corresponding' axisym- 
metric case (Bertola, Rubin & Zeilinger 1989). Similarly, the 
predictions of the triaxial stellar models deviate from those 
in the axisymmetric case: whenever (Tatars > cr ax i sym , then 

^stars ^ ^axisym , 

(ii) Most of the previous modeling assumes that the gas 
is dynamically cold and therefore rotates at the local circu- 
lar speed on the galactic equatorial plane. If in bulges the 
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gas velocity dispersion o" gas is not negligible (e.g. Cinzano 
& van der Marel 1994 hereafter CvdM94; Rix et al. 1995; 
Bertola et al. 1995b), the gas rotates slower than the lo- 
cal circular velocity due to its dynamical pressure support. 
CvdM94 showed explicitly for the E4/S0a galaxy NGC 2974 
that the gas and star kinematics agree taking into account 
for the gas velocity dispersion. Furthermore, if o" gas is com- 
parable to the observed streaming velocity, the spatial gas 
distribution can no longer be modeled as a disc; 

(iii) Forces other than gravity (such as magnetic fields, 
interactions with stellar mass loss envelopes and the hot 
gas component) might act on the ionized gas (e.g. Mathews 
1990). 

In this paper we pursue the second of these explanations 
by building for NGC 4036 dynamical models which take 
into account both for the random motions and the three- 
dimensional spatial distribution of the ionized gas. 

NGC 4036 has been classified S0 3 (8)/Sa in RSA 
(Sandage & Tammann 1981) and S0~ in RC3 (de Vau- 
couleurs et al. 1991). It is a member of the LGG 266 group, 
together with NGC 4041, IC 758, UGC 7009 and UGC 7019 
(Garcia 1993). It forms a wide pair with NGC 4041 with 
a separation of 17' corresponding to 143 kpc at their mean 
redshift distance of 29 Mpc (Sandage & Bedke 1994). In 
The Carnegie Atlas of Galaxies (hereafter CAG) Sandage 
& Bedke (1994) describe it as characterized by an irregu- 
lar pattern of dust lanes threaded through the disc in an 
'embryonic' spiral pattern indicating a mixed SO/Sa form 
(see Panel 60 in CAG). Its total V— band apparent magni- 
tude is Vt — 10.66 mag (RC3). This corresponds to a total 
luminosity Lv = 4.2 ■ 10 10 L0 v at the assumed distance 
of d = 30.2 Mpc. The distance of NGC 4036 was derived 
as d = Vo/Ho from the systemic velocity corrected for the 
motion of the Sun with respect to the centroid of the Local 
Group Vo = 1509±50 km s -1 (RSA) and assuming H = 50 
km s _1 Mpc -1 . At this distance the scale is 146 pc arcsec -1 . 

The total masses of neutral hydrogen and dust in 
NGC 4036 are Mm = 1.7- 10 9 M and M dust = 4.4- 10 5 M 
(Roberts et al. 1991). NGC 4036 is known to have emission 
lines from ionized gas (Bettoni & Buson 1987) and the mass 
of the ionized gas is -Mhii = 7 • 10 4 M (see Sec. [□] for a 
discussion) . 

This paper is organized as follows. In Sec. ^ we 
present the photometrical and spectroscopical observations 
of NGC 4036, the reduction of the data, and the analy- 
sis procedures to measure the surface photometry and the 
major-axis kinematics of stars and ionized gas. In Sec. ^|we 
describe the stellar dynamical model (based on Jeans Equa- 
tions) , and we find the potential due to the stellar bulge and 
disc components starting from the observed surface bright- 
ness of the galaxy. In Sec. ^ we use the derived potential 
to study the dynamics of both the gaseous spheroid and 
disc components, assumed to be composed of collisionless 
cloudlets orbiting as test particles. In Sec. |B| we discuss our 
conclusions. 



2.1.1 Ground-based data 

We obtained an image of NGC 4036 of 300 s in the Johnson 
V— band at the 2.3-m Bok Telescope at Kitt Peak National 
Observatory on December 22, 1995. 

A front illuminated 2048x2048 LICK2 Loral CCD with 
15 x 15 fim 2 pixels was used as detector at the Richtey- 
Chretien focus, f/9. It yielded a flat field of view with a 
10.' 1 diameter. The image scale was Of! 43 pixel" 1 after a 
3x3 pixel binning. The gain and the readout noise were 1.8 
e~ ADU -1 and 8 e~ respectively. 

The data reduction was carried out using standard 
IRAF^| routines. The image was bias subtracted and then 
flat-field corrected. The cosmic rays were identified and re- 
moved. A Gaussian fit to field stars in the resulting image 
yielded a measurement of the seeing point spread function 
(PSF) FWHM of 1'.'7. 

The sky subtraction and elliptical fitting of the galaxy 
isophotes were performed by means of the Astronomical Im- 
ages Analysis Package (AIAP) developed at the Osservato- 
rio Astronomico di Padova (Fasano 1990) . The sky level was 
determined by a polynomial fit to the surface brightness of 
the frame regions not contaminated by the galaxy light, and 
then subtracted out from the total signal. The isophote fit- 
ting was performed masking the bad columns of the frame 
and the bright stars of the field. Particular care was taken 
in masking the dust-affected regions along the major axis 
between 5" and 20". No photometric standard stars were 
observed during the night. For this reason the absolute cali- 
bration was made scaling the total apparent V— band mag- 
nitude to Vt = 10.66 mag (RC3). 

Fig. [I] shows the V— band surface brightness (ftv), ellip- 
ticity (e), major axis position angle (PA), and the cos 49 (124) 
Fourier coefficient of the isophote's deviations from elliptical 
as functions of radius along the major axis. 

For r <; 4" the ellipticity is <; 0.12. Between r ~ 4" 
and r ~ 30" it increases to 0.61. It rises to 0.63 at r ~ 62" 
and then it decreases to 0.56 at the farthest observed radius. 
The position angle ranges from ~ 98° to ~ 67° in the inner 
5". Between r ~ 5" and r ~ 10" it increases to ~ 79° 
and then it remains constant. The cos 4$ coefficient ranges 

between ~ +0.01 and <• 0.02 for r <; 5". Further out 

it peaks to ~ +0.05 at r ~ 18", and then it decreases to 
~ +0.03 for r > 30". The abrupt variation in position angle 
(A PA ~ 30°) observed inside 5" leads to an isophote twist 
that can be interpreted as due to a slight triaxiality of the 
inner regions of the stellar bulge. Anyway this variation has 
to be considered carefully due to the presence of dust pattern 
in these regions which are revealed by HST imaging (see 
Fig-i). 

These results are consistent with previous photometric 
studies of Kent (1984) and Michard (1993) obtained in the 
r— and B— band respectively. Our measurements of ellip- 
ticity follows closely those by Michard (1993). Kent (1984) 
measured the ellipticity and position angle of NGC 4036 
isophotes in r-band for 9" < r < 114". Out to r = 78" 
the r— band ellipticity is lower than our V— band one by 



2 OBSERVATIONS AND DATA ANALYSIS 
2.1 Photometrical observations 



* IRAF is distributed by the National Optical Astronomy Ob- 
servatories which are operated by the Association of Universities 
for Research in Astronomy (AURA) under cooperative agreement 
with the National Science Foundation 
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Figure 1. Ground- based V— band observed surface-brightness, 
cllipticity, position angle and cos 49 coefficient radial profiles of 
NGC 4036. The dashed line in the top panel represents the 
surface-brightness profile of the exponential disc of the best-fit 
kinematical model (fio = 19 .3 mag arcsec -2 , r d = 22'.'0 and 
i = 72°) derived in Sec. 13.2.3 



~ 0.04. The r— band position angle profile differs from our 
only atr~ 16" (e r = 76°) and for r > 78" (e r = 81°). We 
found cos 48 isophote deviation from ellipses to have a radial 
profile in agreement with that by Michard (1993). 



2.1.2 Hubble Space Telescope data 

In addition, we derived the ionized gas distribution in the 
nuclear regions of NGC 4036 by the analysis of two Wide 
Field Planetary Camera 2 (WFPC2) images which were ex- 
tracted from the Hubble Space Telescope archive^]. 

We used a 300 s image obtained on August 08, 1994 
with the F547M filter (principal investigator: Sargent GO- 
05419) and a 700 s image taken on May 15, 1997 with the 
F658N filter (principal investigator: Malkan GO-06785). 

The standard reduction and calibration of the images 
were performed at the STScI using the pipeline- WFPC2 
specific calibration algorithms. Further processing using the 



t Observations with the NASA/ESA Hubble Space Telescope 
were obtained from the data archive at the Space Telescope Sci- 
ence Institute (STScI), operated by AURA under NASA contract 
NAS 5-26555. 



o 

cd 

CO 

a. 



0- 




60 - 



Figure 2. Surface brigthness, ellipticity and position angle radial 
profiles of the nuclear regions of NGC 4036 obtained from a HST 
image (crosses) compared to the ground-based data before (filled 
squ ares) a nd after (open circles) the seeing deconvolution (see 
Sec. 3.2.1 ). The WFPC2 image is in the F547W-band while the 
ground-based ones are in the V— band. 



IRAF STSDAS package involved the cosmic rays removal 
and the alignment of the images (which were taken with 
different position angles). The surface photometry of the 
F547M image was carried out using the STSDAS task EL- 
LIPSE without masking the dust lanes. In Fig. ^|we plot the 
resulting ellipticity (e) and major axis position angle (PA) 
of the isophotes as functions of radius along the major axis. 

The continuum-free image of NGC 4036 (Fig. |^) was 
obtained by subtracting the continuum-band F547M image 
suitably scaled, from the emission-band F658N image The 
mean scale factor for the continuum image was estimated by 
comparing the intensity of a number of 5x5 pixels regions 
near the edges of the frames in the two bandpasses. These 
regions were chosen in the F658N image to be emission free. 
Our continuum-free image reveals that less than the 40% of 
the Ha x [N n] flux of NGC 4036 derives from a clumpy struc- 
ture of about 6" x 2" . The center of this complex filamentary 
structure which is embedded in a smooth emission pattern 
coincides with the position of the maximum intensity of the 
continuum. 



2.2 Spectroscopical observations 

A major-axis (PA= 85°) spectrum of NGC 4036 was ob- 
tained on March 30, 1989 with the Red Channel Spectro- 
graph at the Multiple Mirror Telescope^ as a part of a larger 
sample of 8 SO galaxies (Bertola et al. 1995b). 

The exposure time was 3600 s and the 1200 grooves 



t The MMT is a joint facility of the Smithsonian Institution and 
the University of Arizona. 
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Figure 3. WFPC2 stellar continuum image (upper panel) 
and the continuum-subtracted Ha+[N II] emission image (lower 
panel) of the nucleus of NGC 4036. 



mm -1 grating was used in combination with a l'/25 x 180" 
slit. It yielded a wavelength coverage of 550 A between about 
3650 A and about 4300 A with a reciprocal dispersion of 
54.67 A mm . The spectral range includes stellar absorp- 
tion features, such as the the Ca n H and K lines (AA 3933.7, 
3968.5 A) and the Ca i g-band (A 4226.7 A), and the ion- 
ized gas [O n] emission doublet (AA 3726.2. 3728.9 A). The 
instrumental resolution was derived measuring the a of a 
sample of single emission lines distributed all over the spec- 
tral range of a comparison spectrum after calibration. We 
checked that the measured er's did not depend on wave- 
length, and we found a mean value a = 1.1 A. It corresponds 
to a velocity resolution of ~ 88 km s _1 at 3727 A and ~ 83 
km s _1 at 3975 A. The adopted detector was the 800x800 
Texas Instruments CCD, which has 15x15 /jm 2 pixel size. 
No binning or rebinning was done. Therefore each pixel of 
the frame corresponds to 0.82 A x 0'.'33. 

Some spectra of late-G and early-K giant stars were 
taken with the same instrumental setup for use as velocity 



and velocity dispersion templates in measuring the stellar 
kinematics. Comparison helium-argon lamp exposures were 
taken before and after every object integration. The seeing 
FWHM during the observing night was in between 1" and 
1'.'5. 

The data reduction was carried out with standard pro- 
cedures from the ESO-MIDAS^] package. The spectra were 
bias subtracted, flat-field corrected, cleaned for cosmic rays 
and wavelength calibrated. The sky contribution in the spec- 
tra was determined from the edges of the frames and then 
subtracted. 



2.2. 1 Stellar kinematics 

The stellar kinematics was analyzed with the Fourier Quo- 
tient Method (Sargent et al. 1977) as applied by Bertola et 
al. (1984). The K4III star HR 5201 was taken as template. 
It has a radial velocity of —2.7 km s" 1 (Evans 1967) and 
a rotational velocity of 10 km s _1 (Bernacca & Perinotto 
1970). No attempt was made to produce a master tem- 
plate by combining the spectra of different spectral types, 
as done by Rix & White (1992) and van der Marel et al. 
(1994). The template spectrum was averaged along the spa- 
tial direction to increase the signal-to-noise ratio (S/N). 
The galaxy spectrum was rebinned along the spatial direc- 
tion until a ratio S/N > 10 was achieved at each radius. 
Then spectra of galaxy and template star were rebinned 
to a logarithmic wavelength scale, continuum subtracted 
and endmasked. The least-square fitting of Gaussian broad- 
ened spectrum of the template star to the galaxy spectrum 
was done in the Fourier space over the restricted range of 
wavenumbers [k m i n , k max ] = [5, 200]. In this way we rejected 
the low-frequency trends (corresponding to k < 5) due to the 
residuals of continuum subtraction and the high-frequency 
noise (corresponding to k > 200) due to the instrumental 
resolution. (The wavenumber range is important in particu- 
lar in the Fourier fitting of lines with non-Gaussian profiles, 
see van der Marel & Franx 1993; CvdM94). 

The values obtained for the stellar radial velocity and 
velocity dispersion as a function of radius are given in 



Tab. Al. The table reports the galactocentric distance r 
in arcsec (Col. 1), the heliocentric velocity V (Col. 2) and 
its error 8V (Col. 3) in km s _1 , the velocity dispersion a 
(Col. 4) and its error Sa (Col. 5) in km s _1 . The values for 
the stellar 8V and 5a are the formal errors from the fit in 
the Fourier space. 

The systemic velocity was subtracted from the observed 
heliocentric velocities and the profiles were folded about the 
centre, before plotting. We derive for the systemic heliocen- 
tric velocity a value Vq = 1420 ± 15 km s _1 . Our determina- 
tion is in agreement within the errors with Vq = 1397 ± 27 
km s" 1 (RC3) and V Q = 1382 ± 50 km s" 1 (RSA) derived 
from optical observations too. The resulting rotation curve, 
velocity dispersion profile and rms velocity (%/ v 2 + a 2 ) curve 
for the stellar component of NGC 4036 are shown in Fig. ^. 
The kinematical profiles are symmetric within the error bars 
with respect to the galaxy centre. For r ^ 2" the rotation 
velocity increases almost linearly with radius up to ~ 100 



8 MIDAS is developed and mantained by the European Southern 
Observatory 
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Figure 4. Stellar kinematics along the major axis of NGC 4036. 
The folded rotation velocity curve (top panel), velocity dispersion 
profile (middle panel) and rms velocity curve (bottom panel) are 
shown in km s — 1 . Open squares and filled triangles represent data 
derived for the approaching W and receding E side respectively. 



km s , remaining approximatively constant between 2" 
and 4". Outwards it rises to the farthest observed radius. 
It is ~ 180 km s _1 at 9", ~ 220 km s _1 at 21" and ~ 260 
km s _1 at 29". The velocity dispersion a ~ 210 km s _1 in 
the centre and at r ~ 4" with a 'local minimum of ~ 180 at 
r = 2". Further out it declines to values ^ 120 km s . 



2.2.2 Ionized gas kinematics 

To determine the ionized gas kinematics we studied the [O n] 
(AA 3726.2, 3728.9 A) emission doublet. In our spectrum the 
two lines are not resolved at any radius. We obtained smooth 
fits to the [O n] doublet using a two-steps procedure. In the 
first step the emission doublet was analyzed by fitting a dou- 
ble Gaussian to its line profile, fixing the ratio between the 
wavelengths of the two lines and assuming that both lines 
have the same dispersion. The intensity ratio of the two lines 
depends on the state (i.e. electron density and temperature) 
of the gas (e.g. Osterbrock 1989). We found a mean value of 
[O n] A3726.2/[0 n] A3728.9 = 0.8 ±0.1 without any signi- 
ficative dependence on radius. The electron density derived 




10 20 
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Figure 5. Ionized gas kinematics along the major axis of 
NGC 4036. The folded rotation velocity curve (top panel), ve- 
locity dispersion profile (middle panel) and rms velocity curve 
(bottom panel) are shown in km s . Open squares and filled tri- 
angles represent data derived for the approaching W and receding 
E side respectively. The velocity error bars indicate the formal er- 
ror of double Gaussian fit to the [O II] doublet. The error bars in 
the velocity dispersions take also into account for the subtraction 
of the instrumental dispersion. 



from the obtained intensity ratio of [O n] lines is in agree- 
ment with that derived (at any assumed electron temper- 
ature, see Osterbrock 1989) from the intensity ratio of the 
[S n] lines ([S n] A6716.5/[S n] A6730.9=1.23) found by Ho, 
Filippenko & Sargent (1997). In the second step we fitted 
the line profile of the emission doublet by fixing the inten- 
sity ratio of its two lines at the value above. At each radius 
we derived the position, the dispersion and the uncalibrated 
intensity of each [O n] emission line and their formal errors 
from the best-fitting double Gaussian to the doublet plus a 
polynomial to its surrounding continuum. The wavelength 
of the lines' centre was converted into the radial velocity 
and then the heliocentric correction was applied. The lines' 
dispersion was corrected for the instrumental dispersion and 
then converted into the velocity dispersion. 

The measured kinematics for the gaseous component 
in NGC 4036 is given in Tab. |a|. The table contains the 
galactocentric distance r in arcsec (Col. 1), the heliocentric 
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Figure 6. The intensity of the [O II] A3726.2 line as function of 
radius in NGC 4036. Since the spectrum was not flux calibrated 
the scale has an arbitrary zero point. The thick continuous lines 
connect the measurements along the W and E side respectively. 
The vertical bar on the left-hand side of the panel represents the 
typical error for the data. The dashed and dotted lines represent 
respectively the profile of the gas eou s i? 1 / 4 spheroid and gaseous 
exponential disc derived in Sec. [l. l| . The thin continuous line is 
the fit to the data. 



velocity V (Col. 2) and its error 5V (Col. 3) in km s _1 the 
velocity dispersion a (Col. 4) and its errors 5a+ and 5a- 
(Cols. 5 and 6) in km s _1 . The gas velocity errors 8V are 
the formal errors for the double Gaussian fit to the [O n] 
doublet. The gas velocity dispersion errors 8a + and <5<r_ 
take also account for the subtraction of the instrumental 
dispersion. 

The rotation curve, velocity dispersion profile and rms 
velocity curve for the ionized gas component of NGC 4036 
resulting after folding about the centre are shown Fig. || 
The [O n] A3726.2 intensity profile as a function of radius is 
plotted in Fig. [| The gas rotation tracks the stellar rotation 
remarkably well. They are consistent within the errors to one 
another. The gas velocity dispersion has central dip of ~ 160 
km s _1 with a maximum of ~ 220 km s~ x at r ~ 2". It re- 
mains higher than 100 km s" 1 up to ~ 4" before decreasing 
to lower values. The velocity dispersion profile appears to 
be less symmetric than the rotation curve. Indeed between 
2" and 5" the velocity dispersion measured onto the E side 
rapidly drops from its observed maximum to ~ 50 km s _1 , 
while in the W side it smoothly declines to ~ 140 km s _1 . 
Errors on the gas velocity dispersion increase at large radii 
as the gas velocity dispersion becomes comparable to the 
instrumental dispersion. 



2.2.3 Comparison with kinematical data by Fisher (1997) 

The major-axis kinematics for stars and gas we derived for 
NGC 4036 are consistent within the errors with the mea- 



surements of Fisher (1997, hereafter F97). The only excep- 
tion is represented by the differences of 20%-30% between 
our and F97 stellar velocity dispersions in the central 8". In 
these regions F97 finds a flat velocity dispersion profile with 
a plateau at a ~ 170 km s _1 . To measure stellar kinemat- 
ics he adopted the Fourier Fitting Method (van der Marel & 
Franx 1993) directly on the line-of-sight velocity distribution 
derived with Unresolved Gaussian Decomposition Method 
(Kuijken & Merrifield 1993). For \r\ <, 8" the NGC 4036 
line profiles are asymmetric (displaying a tail opposite to 
the direction of rotation) and flat-toped, as result from the 
/13 and hi radial profiles. The /13 and /14 parameters mea- 
sure respectively the asymmetric and symmetric deviations 
of the line profile from a Gaussian (van der Marel & Franx 
1993; Gerhard 1993). For NGC 4036 the h 3 term anticorre- 
lates with v, rising to ~ +0.1 in the approaching side and 
falling to ~ —0.1 in the receding side. The /14 term exhibits 
a negative value (~ —0.03). 



3 MODELING THE STELLAR KINEMATICS 
3.1 Modeling technique 

We built an axisymmetric bulge-disc dynamical model for 
NGC 4036 applying the Jeans modeling technique intro- 
duced by Binney, Davies & Illingworth (1990), developed 
by van der Marel, Binney & Davies (1990) and van der 
Marel (1991), and extended to two-component galaxies by 
CvdM94 and to galaxies with a DM halo by Cinzano (1995) 
and Corsini et al. (1998). For details the reader is referred 
the above references. 

The main steps of the adopted modeling are (i) the 
calculation of the bulge and disc contribution to the poten- 
tial from the observed surface brightness of NGC 4036; (ii) 
the solution of the Jeans Equations to obtain separately the 
bulge and disc dynamics in the total potential; and (iii) the 
projection of the derived dynamical quantities onto the sky- 
plane taking into account seeing effects, instrumental set-up 
and reduction technique to compare the model predictions 
with the measured stellar kinematics. In the following each 
invidual step is briefly discussed: 

(i) We model NGC 4036 with an infinitesimally thin ex- 
ponential disc in its equatorial plane. The disc surface mass 
density is specified for any inclination i, central surface 
brightness po, scale length rd and constant mass-to- light ra- 
tio (M/L)d- The disc potential is calculated from the sur- 
face mass density as in Binney & Tremaine (1987). The 
limited extension of our kinematical data (measured out 

r < 30" ~ 0.5 .RopiPj]) prevents us to disentangle in the 
assumed constant mass-to-light ratios the possible contribu- 
tion of a dark matter halo. 

The surface brightness of the bulge is obtained by subtract- 
ing the disc contribution from the total observed surface 
brightness. The three-dimensional luminosity density of the 
bulge is obtained deprojecting its surface brightness with an 
iterative method based on the Richardson-Lucy algorithm 
(Richardson 1972; Lucy 1974). Its three-dimensional mass 

^ The optical radius R pt is the radius encompassing the 83% of 
the total integrated light. 
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density is derived by assuming a constant mass-to-light ra- 
tio (M/L)b. The potential of the bulge is derived solving the 
Poisson Equation by a multipole expansion (e.g. Binney & 
Tremaine 1987). 

(ii) The bulge and disc dynamics are derived by sepa- 
rately solving the Jeans Equations for each component in 
the total potential of the galaxy. For both components we 
assume a two integral distribution function of the form 
/ = f(E,L z ). It implies that the vertical velocity disper- 
sion <r 2 is equal to second radial velocity moment o\ and 
that vrv z — 0. Therefore the Jeans Equations becomes a 
closed set, which can be solved for the unknowns and 
a 2 R — <r 2 . Other assumptions to close the Jeans Equations 
are also possible (e.g. van der Marel & Cinzano 1992). 

For the bulge we made the same hypotheses of Binney et 
al. (1990). A portion of the second velocity moment w? is 
assigned to bulge streaming velocity v$ following Satoh's 
(1980) prescription. 

For the disc we made the same hypotheses of Rix & White 
(1992) and CvdM94. The second radial velocity moment o\ 
in the disc is assumed to fall off exponentially with a scale 
length R a from a central value of aj . The azimuthal ve- 
locity dispersion in the disc is assumed to be related 
to a\ according to the relation from epicyclic theory [cfr. 
Eq. (3-76) of Binney & Tremaine (1987)]. As pointed out 
by CvdM94 this relation may introduce systematic errors 
(Kuijken & Tremaine 1992; Evans & Collett 1993; Cudde- 
ford & Binney 1994). The disc streaming velocity v$ (i.e. 
the circular velocity corrected for the asymmetric drift) is 
determined by the Jeans equation for radial equilibrium. 

(iii) We projected back onto the plane of the sky (at the 
given inclination angle) the dynamical quantities of both the 
bulge and the disc to find the line-of-sight projected stream- 
ing velocity and velocity dispersion. We assumed that both 
the bulge and disc have a Gaussian line profile. At each ra- 
dius their sum (normalized to the relative surface brightness 
of the two components) represents the model-predicted line 
profile. As in CvdM94 the predicted line profiles were con- 
volved with the seeing PSF of the spectroscopic observations 
and sampled over the slit width and pixel size to mimic the 
observational spectroscopic setup. We mimicked the Fourier 
Quotient method for measuring the stellar kinematics fitting 
the predicted line profiles with a Gaussian in the Fourier 
space to derive the line-of-sight velocities and velocity dis- 
persions for the comparison with the observed kinematics. 
The problems of comparing the true velocity moments with 
the Fourier Quotient results were discussed by van der Marel 
& Franx (1993) and by CvdM94. 



3.2 Results for the stellar component 

3.2.1 Seeing- deconvolution 



The modeling technique described in Sec. 3.1 derives 



the three-dimensional mass distribution from the three- 
dimensional luminosity distribution inferred from the ob- 
served surface photometry by a fine-tuning of the disc pa- 
rameters leading to the best fit on the kinematical data. We 
performed a seeing-deconvolution of the V— band image of 
NGC 4036 to take into account the seeing effects on the mea- 
sured photometrical quantities (surface-brightness, elliptic- 
ity and cos 46 deviation profiles) used in the deprojection of 



the two-dimensional luminosity distribution. We obtained 
a restored NGC 4036 image through an iterative method 
based on the Richardson-Lucy algorithm (Richardson 1972; 
Lucy 1974) available in the IRAF package STSDAS. We 
assumed the seeing PSF to be a circular Gaussian with a 
FWHM= l'.'7. Noise amplification represents a main back- 
draw in all Richardson-Lucy iterative algorithms and the 
number of iterations needed to get a good image restoration 
depends on the steepness of the surface-brightness profile 
(e.g. White 1994). After 6 iterations, we did not notice any 
more substantial change in the NGC 4036 surface-brightness 
profile while the image became too noisy. Therefore we de- 
cided to stop at the sixth iteration. The surface-brightness, 
ellipticity, position angle radial profiles of NGC 4036 after 
the seeing-deconvolution are displayed in Fig. [] compared 
with the HST and the unconvolved ground-based photome- 
try. The raises found at small radii for the surface brightness 
(r ~ 0.2 mag arcsec -2 ) and for the ellipticity (r ~ 0.05) are 
comparable in size to the values found by Peletier et al. 
(1990) for seeing effects in the centre of ellipticals. 



3.2.2 Bulge-disc decomposition 

We performed a standard bulge-disc decomposition with a 
parametric fit (e.g. Kent 1985) in order to find a starting 
guess for the exponential disc parameters to be used in 
the kinematical fit. We decomposed the seeing-deconvolved 
surface-brightness profile on both the major and the minor 
axis as the sum of an 7? 1//4 bulge, having surface-brightness 
profile 



MO = fJ. e + 8.3268 



(1) 



plus an exponential disc, having surface-brightness profile 



MO = Mo + 1-0857 



Td 



(2) 



We assumed that the minor-axis profile of each component 
is the same as the major-axis profile, but scaled by a factor 
1 — e = b/a. A least-squares fit of the photometric data 
provides fj, e , r e and e of the bulge, /Lto, r<j of the disc, and 
the galaxy inclination i (see Tab. hi for the results). 

Kent (1984) measured surface photometry of NGC 4036 
in the r— band. He decomposed the major- and minor-axis 
profiles in an R 1 ^ 4 bulge and in an exponential disc (Kent 
1985). A rough but useful comparison between the bulge 
and disc parameters resulting from the two photometric 
decompositions is possible by a transformation from r- to 
V— band of Kent's data. We derived from Kent's surface 
brightness profile along the major axis of NGC 4036 its (ex- 
trapolated) total magnitude vt = 10.56L0.02 corresponding 
to Vt — Tt = 0.13 ± 0.11. Then we converted the fio and /j, e 
values from the r- to the V— band (see Tab. [I]). The dif- 
ferences between the best-fit parameters obtained from the 
two decompositions are lower than 10% except for bulge el- 
lipticity (our value is ~ 64% than the Kent's one). These 
discrepancies are consistent with the differences in the slope 
of the two surface-brightness profiles, with differences be- 
tween the Kent's approach to correct for the seeing effects 
[Kent (1985) convolved the theoretical bulge and disc pro- 
files with the observed Gaussian seeing profile] and with the 
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Table 1. The bulge-disc decomposition parameters 
bulge disc 



this paper 
Kent (1985) 



/tie 

20.4 
20.7 



12"8 
13'.' 1 



0.11 
0.07 



Mo 
18.7 
18.6 



rd i 
22'.' 1 74? 9 
21'.' 71° 9 



Note. fi e and no are given in V— mag arcsec 2 . Kent's (1984) sur- 
face brightnesses have been converted from the r- to the V— band. 

uncertainties of the conversion from r- to V— band of Kent's 
surface brightnesses. 



3.2.3 Modeling results 

We looked for the disc parameters leading to the best fit 
of the observed stellar kinematics, using as starting guess 
those resulting from our best-fit photometric decomposition, 
with |/io — 1| < 18.7 mag arcsec -2 , \rd — 5" < 22'.' 1, and 
\i — 5° | < 74? 9. For any exponential disc of fixed /io, rd, i (in 
the investigated range of values), we subtracted its surface 
brightness from the total seeing-deconvolved surface bright- 
ness. The residuum surface brightness was considered to be 
contributed by the bulge. Being /io, rd and i correlated, a 
bulge component with a surface-brightness profile consis- 
tent with that resulting from the photometric decomposi- 
tion, was obtained only by taking exponential discs charac- 
terized by large /io values in combination with large values 
of rd (i.e. larger but fainter discs) or by small /io and small 
rd (i.e. smaller but brighter discs). After subtracting the 
disc contribution to the total surface brightness, the three- 
dimensional luminosity density of the bulge has been ob- 
tained after seven Richardson-Lucy iterations starting from 
a fit to the actual bulge surface brightness with the projec- 
tion of a flattened Jaffe profile (1983). The residual surface 
brightness (A/i = /imodel — Mobs) after each iteration and 
the final three-dimensional luminosity density profiles of the 
spheroidal component of NGC 4036 along the major, the mi- 
nor and two intermediate axis is plotted for the kinematical 
best-fit model in Fig. ^ up to 100" (the stellar and ionized- 
gas kinematics are measured to ~ 30"). 



We applied the modeling technique described in Sec. 3.1 
considering only models in which the bulge is an oblate 
isotropic rotator (i.e. k = 1) and in which the bulge and the 
disc have the same mass-to-light ratio (M/L)b = (M/L)d- 
The M/L determines the velocity normalization and was 
chosen (for each combination of disc parameters) to opti- 
mize the fit to the rms velocity profile. The best-fit model to 
the observed major-axis stellar kinematics is obtained with 
a mass-to-light ratio M/Lv= 3.42 Mq/Lqv and with an 
exponential disc having a surface-brightness profile 



Vu(r) = 19.3 + 1.0857 (^) 



mag arcsec 



(3) 



(represented by the dashed line in the upper panel of Fig. |l|), 
a radial velocity dispersion profile 



crfl(r) = 155 e r ' 27 ' 4 km s 1 , 



(4) 



(where the galactocentric distance r is expressed in arcsec) 
and an inclination i = 72°. Fig. ^ shows the comparison 
between the rotation curve, the dispersion velocity profile 
and the rms velocity curve predicted by the best-fit model 





3 10 30 100 
R/ arcsec 



Figure 7. The deprojection of the surface brightness of the 
spheroidal component of NGC 4036 for the best-fit kinematical 
model. The right panel shows the residual surface brightness A/i 
in mag arcsec - 2 after each Lucy iteration (dashed lines) from 
an initial Jaffe fit to the actual NGC 4036 bulge brightness. The 
residuals are shown for four axes (major through minor axis: top 
through bottom). For each set of curves: the solid line corresponds 
to the projected adopted model for the 3-D luminosity density and 
the dotted line corresponds to a perfect deprojection. At each in- 
teraction if the model is brighter than the galaxy, the correspond- 
ing dashed or continuous curve is below the dotted line. In the 
left panel the final three-dimensional luminosity density profile 
of the spheroidal component of NGC 4036 right panel is given 
in units of 1O 1O L0 kpc - 3 for the same four axes (minor through 
major axis: innermost through outermost curve). 



(solid lines) with the observed stellar kinematics along the 
major axis of NGC 4036. The agreement is good. 

The derived V— band luminosities for bulge and disc 
are L b = 2.8 • 10 10 L ev and L d = 1.4 • 10 10 L ev . They 
correspond to the masses M b = 9.8 • 10 10 M Q and M d = 
4.8 • 10 10 Mq. The total mass of the galaxy (bulge+disc) is 
M T = 14.5 • 10 10 M Q . The disc-to-bulge and disc-to-total 
V— band luminosity ratios are Lt/Ld = 0.58 and Ld/Lr = 
0.36. The disc-to-total luminosity ratio as function of the 
galactocentric distance is plotted in Fig. ^. 



3.2.4 Uncertainty ranges for the disc parameters 

The uncertainty ranges for the disc parameters are 19.1 
/io ^ 19.5 mag arcsec" 2 , 20" <, r d & 24" and 71° <, i <, 73°. 
In Fig. [h^ the continuous and the dotted line show the kine- 
matical profiles predicted for two discs with the same in- 
clination (i = 72°) and the same total luminosity of the 
best-fit disc but with a smaller (/io = 19.1 mag arcsec -2 , 
rd — 20'.'0) or greater (jjlq — 19.5 mag arcsec -2 , ra = 24'.' 0) 
scale length respectively. In Fig. jH] the continuous and the 
dotted lines correspond to two discs with the same scale 
length (rd = 22'.' 0) of the best-fit disc but with a lower 
(no = 19.2 mag arcsec -2 , i = 71°) or a higher (fj,o = 19.4 
mag arcsec -2 , i = 73°) inclination. Their total luminosity 
is respectively ~ 14% lower and ~ 16% higher then that of 
the best-fit disc. The fit of these models to the data are also 
acceptable so we estimate a ~ 15% error in the determina- 
tion of the disc luminosity and mass. The good agreement 
with observations obtained with k = 1, (M/L)b = (M/L)d, 
and without a dark matter halo causes us to not investigate 
models with k 7^ 1, (M/L)b 7^ (M/L)d or with radially in- 
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Figure 8. Comparison of the model predictions to the stellar 
major-axis kinematic for NGC 4036. Data points are as in Fig. Q 
The solid curves represent the velocity (upper panel), velocity 
dispersion (middle panel) and rms velocity (lower panel) radial 
profiles of the best-fit model (in which the disc parameters are 
fig = 19.3 mag arcsec -2 , r d = 22'.'0 and i = 72°). 



creasing mass-to-light ratios. Therefore we cannot exclude 
them at all. 

The V— band luminosity (after inclination correction) 
of the exponential disc corresponding to the model best-fit 
to the observed stellar kinematics is ~ 75% than that of 
the disc obtained from the parametric fit of the surface- 
brightness profiles. The differences between the disc pa- 
rameters derived from NGC 4036 photometry (/xo = 18.7 
mag arcsec -2 , rd = 22 / .'l, i = 74?9) and from the stellar 
kinematics (no = 19.3 mag arcsec -2 , = 22'.'0, i = 72°) 
are expected. In fitting the NGC 4036 surface-brightness 
profiles the bulge is assumed (i) to be axisymmetric; (ii) 
to have a R 1 ^ 4 profile; (iii) and constant axial ratio (i.e. 
its isodensity luminosity spheroids are similar concentric el- 
lipsoids). We adopted this kind of representation for the 
bulge component only to find rough bounds on the expo- 
nential disc parameters to be used in the kinematical fit. 
Often bulges have neither an _R 1//4 law profile (e.g. Burstein 
1979; Simien & Michard 1984) neither perfectly elliptical 
isophotes (e.g. Scorza & Bender 1990). Therefore in model- 
ing the stellar kinematics the three-dimensional luminosity 
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Figure 9. Radial profile of the fraction of the light contributed 
by the stellar disc along the major axis of NGC 4036 as obtained 
by the best-fit model. The distance from the center is expressed 
in arcsec. 
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Figure 10. As in Fig 



mag arcsec -2 , r d = 20'.'0 and i = 
for a disc with fio = 19.5 mag arcsec 
(dotted line). 



but for a disc with fio = 19.1 
72° (continuous line) and 
" 2 , r d = 24'.'0 and i = 72° 



density jt of the stellar bulge is assumed (i) to be oblate 
axisymmetric; but (ii) not to be parametrized by any an- 
alytical expression; (iii) nor to have isodensity luminosity 
spheroids with constant axial ratio. The flattening of the 
spheroids increases with the galactocentric distance, as it 
appears from the three-dimensional luminosity density pro- 
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Figure 11. As in Fig. ^ but for a disc with fio = 19.2 
mag arcsec -2 , = 22'.'0 and i = 71° (continuous line) and 
for a disc with /iq = 19 A mag arcsec -2 , r d = 22" and i = 73° 
(dotted line). 



files plotted in Fig. |t] along different axes onto the meridional 
plane of NGC 4036. This flattening produces an increasing 
of the streaming motions for the bulge component assumed 
to be an isotropic rotator. 



3.2.5 Modeling results with stellar kinematics by Fisher 
(1997) 

The bulge and disc parameters found for our best-fit model 
reproduce also the major-axis stellar kinematics by F97. 
Rather than choosing an 'ad hoc' wavenumber range to em- 
ulate F97 analysis technique the fit to the predicted line 
profile was done in ordinary space and not in the Fourier 
space as done to reproduce our Fourier Quotient measure- 
ments. The good agreement of the resulting kinematical pro- 
files with F97 measurements are shown in Fig. gjl 



MODELING OF THE IONIZED GAS 
KINEMATICS 
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Figure 12. Comparison between the predictions of the same 
model of Fig. pi with the stellar major-axis kinematics for NGC 
4036 by Fisher (1997). Data points are as in Fig. ^ The solid 
curves represent the velocity v, velocity dispersion <r, /13 and /14 
radial profiles of the best-fit model taking into account for Fisher's 
data reduction technique. 



4.1 Modeling technique 

At small radii both the ionized gas velocity and velocity dis- 
persion are comparable to the stellar values, for r < 9" and 
r < 5" respectively. This prevented us to model the ionized 
gas kinematics by assuming the gas had settled into a disc 
component as done by CvdM94 for NGC 2974. Moreover 
a change in the slope of the [O 11] A3726.2 intensity radial 
profile (Fig. ^) is observed at r ^ 8", its gradient appears to 
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be somewhat steeper towards the centre. The velocity dis- 
persion and intensity profiles of the ionized gas suggest that 
it is distributed into two components (see also the distinct 
central structure in HST Hq+[N ii] image, Fig. ^): a small 
inner spheroidal component and a disc. 

We built a dynamical model for the ionized gas with a 
dynamically hot spheroidal and in a dynamically colder disc 
component. 

The total mass of the ionized gas Mhii is negligible 
and the total potential is set only by the stellar component. 
The mass of the ionized gas can be derived from optical 
recombination line theory (see Osterbrock 1989) by the Ha 
luminosity. For a given electron temperature T c and density 
N e , the H ii mass is given by 

Mhii = (Ln a m H /N c ) / (4Trj Ha /N c N p ) (5) 

where Lua is the Ha luminosity, ran is the mass of the hy- 
drogen atom, jHa is the Ha emissivity, and A^p are the proton 
density (Tohline & Osterbrock 1976). The Ha luminosity of 
NGC 4036 scaled to the adopted distance is Lho = 5.6- 10 39 
Lq (Ho et al. 1997). The term 4-Kju a /N e Np is insensitive to 
changes of N e over the range 10 2 - 10 6 cm . It decreases 
by a factor 3 for changes of T c over the range 5 ■ 10 3 °K 
- 2 ■ 10 4 °K (Osterbrock 1989). For an assumed tempera- 
ture T c = 10 4 °K, the electron density is estimated to be 
N c — 2 ■ 10 2 cm -3 from the [S n] ratio found by Ho et al. 
(1997) implying M H n = 7- 10 4 M . 

For the gaseous spheroid and disc we made two different 
sets of assumptions based on two different physical scenarios 
for the gas cloudlets. 

4-1.1 Long-lived gas cloudlets (model A) 

In a first set of models we described the gaseous component 
by a set of collisionless cloudlets in hydrostatic equilibrium. 
The small gaseous 'spheroid' is characterized by a density 
distribution and flattening different from those of stars. Its 
major-axis luminosity profile was assumed to follow an 7? 1//4 
law. Adopting for it the Ryden's (1992) analytical approxi- 
mation, we obtained the three-dimensional luminosity den- 
sity of the gaseous spheroid. The flattening of the spheroids 
q was kept as free parameter. To derive the kinematics of the 
gaseous spheroid we solved the Jeans Equations under the 
same assumptions made in Sec. ^| for the stellar spheroid. 
In particular the streaming velocity v$ of gaseous bulge is 
derived from the second azimuthal velocity moment us- 
ing Satoh's (1980) relation. For the ionized gaseous disc we 
solved the Jeans Equations under similar assumptions made 
in Sec. ^| for the stellar disc. Specifically, we assumed that 
the gaseous disc (i) has an exponential luminosity profile; 
(ii) is infinitesimally thin; (iii) has an exponentially decreas- 
ing a R ; (iv) has a z = a R ; (v) has satisfying the epicyclic 
relation. 

4-1.2 Gas cloudlets 'just' shed by the stars (model B) 

In a second set of models we assumed that the emission ob- 
served in the gaseous spheroid and disc arise from material 
that was recently shed from stars. Different authors (Bertola 
et al. 1984; Fillmore et al. 1986; Kormendy & Westpfahl 
1989; Mathews 1990) suggested that the gas lost (e.g. in 



planetary nebulae) by stars was heated by shocks to the 
virial temperature of the galaxy within 10 4 years, a time 
shorter than the typical dynamical time of the galaxy. Hence 
in this picture the ionized gas and the stars have the same 
true kinematics, while their observed kinematics are differ- 
ent due to the line-of-sight integration of their different spa- 
tial distribution. Differences between the radial profiles of 
the gas emissivity and the stellar luminosity may be ex- 
plained if both the gas emission process and the efficiency of 
the thermalization process show a variation with the galac- 
tocentric distance. The three-dimensional luminosity density 
of the spheroid is derived as in model A and the luminosity 
density profile of the disc is assumed to be exponential. 

In both cases, the kinematics of the gaseous spheroid 
and disc were projected on the sky-plane to be compared to 
the observed ionized gas kinematics. Assuming a Gaussian 
line profile for both the gaseous components we derived the 
total line profile (which depends on the relative flux of the 
two components). As for the stellar components, we con- 
volved the line profiles obtained for the ionized gas with the 
seeing PSF and we sampled them over the slit-width and 
pixel size to mimic the observational setup. This procedure 
is particularly important for the modeling of the observed 
kinematics near the centre. As a last step (in mimicking the 
measuring technique of the gaseous kinematics) we fitted a 
single Gaussian to the resulting line profile, after taking into 
account for the instrumental line profile. 



4.2 Results for the gaseous component 

We decomposed the [O n] A3726.2 intensity profile as the 
sum of an R 1 ^ 4 gaseous spheroid and an exponential gaseous 
disc. A least-squares fit of the observed data was done for 
r > 3" to deal with seeing effects (Fig. ^|). The gas spheroid 
resulted to be the dominating component up to r ~ 8" be- 
yond the bright emission in Fig. |^. We derived the effective 
radius of the gaseous spheroid r e , gas = 0'.'5 ± O'.'l, the scale 
length of the gaseous disc r^gas = 29'.'8 ± (f.'9, and the ra- 
tio between the effective intensity of the spheroid and the 
central intensity of the disc I £> gas //o, gas = 718l*5 3 . The un- 
certainties on the resulting parameters have been estimated 
by a separate decomposition on each side of the galaxy. With 
these parameters we app lied the models for the gas kinemat- 
ics described in Sec. 4.1 

Since in both models A and B the stellar density ra- 
dial profile differs from the gas emissivity radial profile, 
it is interesting to check what is the relation between the 
three-dimensional stellar density p s tar(^i) and the three- 
dimensional gas emissivity f gas (7?). After deprojection we 
find they are related by the following relation 

(R) OC Pstar (R) (6) 

in the range of galactocentric distances between r ~ 2" and 
r ~ 10". In a fully ionized gas the recombination rate is pro- 
portional to the square of the gas density (e.g., Osterbrock 
1989), this is a power-law relation quite similar to Eq. ^. 

For models A and B the best-fit to the observed gas 
kinematics are in both cases obtained with a spheroid flat- 
tening q — 0.8, and are plotted respectively in Fig. [l3] and 
Fig. 0. 

The best result is obtained with model B. A simple estimate 
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Figure 13. Comparison of the predictions of model A (dashed 
curve) to the ionized gas kinematics observed along the major 
axis of NGC 4036. The data points are as in Fig. g. The dotted 
curves represent the seeing-convolved circular velocity curve and 
zero velocity-dispersion profile in the galaxy meridional plane. 



of the errors on the model due to the uncertainties of bulge- 
disc decomposition of the radial profile of [O n] A3726.2 
emission line can be inferred by comparing the model pre- 
dictions based on the separate decomposition on each side 
of the galaxy. We find a maximum difference of 5% for 
4" < r < 10" between the gas velocities and velocity dis- 
persion predicted using the two different bulge-disc decom- 
positions of the [O n] A3726.2 intensity profile. For model 
A the assumption of Satoh's (1980) relation fails in repro- 
ducing the observed gas kinematics for r ^ 6", where the 
emission lines intensity profile is dominated by the gaseous 
spheroidal component. However the R 1 ' 4 extrapolation of 
the density profile in the inner 3" overestimates the density 
gradient in this region (as it appears also from the HST im- 
age, Fig. ^) which could produce an exceeding asymmetric 
drift correction. 



4-2.1 Modeling results with gas kinematics by Fisher 
(1997) 

We also applied our models A and B to the ionized gas 
kinematics and to the [O in] A5006.9 intensity radial profile 
measured by F97 along the major axis of NGC 4036. The 
best-fit to F97 data for models A and B are obtained with a 
spheroid flattening q = 0.8. They are shown in Fig. |l5| and 
Fig. |l6| respectively. 

In this case the dynamical predictions of model B (even 
if they give better results than those of model A) are not 
able to reproduce the F97's kinematics in the radial range 
between 3" and 10". The differences between the predicted 
and the measured kinematics rise up to 80 km s _1 in ve- 
locities and to 40 km s _1 in velocity dispersions at r ~ 6". 
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Figure 14. As in Fig. ^ but for model B. 
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Figure 15. As in Fig. |l3| but for the gas kinematics measured 
by Fisher (1997). Open triangles and filled squares represent data 
derived for the approaching W and receding E side respectively. 



Nevertheless for 0" < r < 10" this model agrees better with 
the observations than if the gas is assumed to be on circular 
orbits (see the dotted lines in Fig. |lri| ). If this is the case 
the maximum differences with the measured kinematics are 
large as 130 km s _1 in velocity and as 110 km s _1 in velocity 
dispersion for r ~ 4". 
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Figure 16. As in Fig. |l^ but for model B. 



4.3 Do drag forces affect the kinematics of the 
gaseous cloudlets? 

Considerable differences exist between the ionized gas kine- 
matics recently measured by F97 and the velocity and ve- 
locity dispersion profiles predicted even by the best model 
in the bulge-dominated region between r ~ 4" and r ~ 10". 
This suggests that other phenomena play a role in deter- 
mining the dynamics of the spheroid gas (see item (iii) in 
Sect. hi). The discrepancy between model and observations 
could be explained by accounting for the drag interaction 
between the ionized gas and the hot component of the in- 
terstellar medium. 

In the scenario of evolution for stellar ejecta in elliptical 
galaxies outlined by Mathews (1990), a portion of the gas 
shed by stars (e.g., as stellar winds or planetary nebulae) 
undergoes to an orbital separation from its parent stars by 
the interaction with ambient gas, after an expansion phase 
and the attainment of pressure equilibrium with the envi- 
ronment, and before its disruption by various instabilities. 
Indeed to explain the luminosity of the optical emission lines 
measured for nearby ellipticals, Mathews (1990) estimated 
that the ionized gas ejected from the orbiting stars merges 
with the hot interstellar medium in at least tu{ e ~ 10 6 yr. 
If the gas clouds have at the beginning the same kinematics 
of their parent stars, this lifetime is sufficiently long to let 
the gaseous cfouds (which start with the same kinematics of 
their parent stars) to acquire an own kinematical behaviour 
due to the deceleration produced by the drag force of the 
intersteilar diffuse medium. The lifetime of the ionized gas 
nebulae is shorter than 10 4 — 10 5 years if magnetics effects 
on gas kinematics are ignored (as in our model B). 

To have some qualitative insights in understanding the 
effects of a drag force on the gas kinematics we studied the 
case of a gaseous nebula moving in the spherical potential 



generated by an homogeneous mass distribution of density 
p and which, starting onto a circular orbit, is decelerated by 
a drag force 



A- 



drag — 



„ 2 V 
V — 

m v 



(8) 



where m and v are the mass and the velocity of the gaseous 
cloud. Following Mathews (1990), the constant kd rag is given 
by 

^ 3 n m 

ft drag ~ "7 



0) 



where n is the density of the interstellar medium, n eq and 
a eq are respectively the density and the radius of the gaseous 
nebula when the equilibrium is reached between the inter- 
nal pressure of the cloud and the external pressure of the 
interstellar medium. The ratio n/n eq ~ 10~ 3 at any galac- 
tic radius and therefore the ratio kdrag /m depends on the 
nebula radius a eq (Mathews 1990). 

The equations of motion of the nebula expressed in 
plane polar coordinates (r, tjj) in which the centre of attrac- 
tion is at r = and ip is the azimuthal angle in the orbital 
plane are 



■r^ 2 = 4 T Gpr + ^f 2 (r>0) (10) 
3 m 



ri: + 2iPr = -^^-rip 2 (r>0) (11) 
m 

We numerically solved the Eqs. [l(j and jll] with the 
Runge-Kutta method (Press et al. 1986) to study the time- 
dependence of the radial and tangential velocity components 
r and rijj of the nebula. We fixed the potential assuming 
a circular velocity of 250 km s _1 at r = 1 kpc. Follow- 
ing Mathews (1990) we took an equilibrium radius for the 
gaseous nebula a eq — 0.37 pc. The results obtained for dif- 
ferent times in which the drag force decelerate the gaseous 
clouds are shown in Fig. [jj]. 

It results that tjj < and r > 0: the clouds spiralize 
towards the galaxy centre as expected. Moreover the drag 
effects are greater on faster starting clouds and therefore 
negligible for the slowly moving clouds in the very inner 
region of NGC 4036. 

If the nebulae are homogeneously distributed in the gaseous 
spheroid only the tangential component rip of their velocity 
contribute to the observed velocity. No contribution derives 
from the radial component r of their velocities. In fact for 
each nebula moving towards the galaxy centre which is also 
approaching to us we expect to find along the line-of-sight 
a receding nebula which is falling to the centre from the 
same galactocentric distance with an opposite line-of-sight 
component of its r. 

The radial components of the cloudlet velocities (typ- 
ically of 30-40 km s~\ see Fig. |l7| ) are crucial to explain 
the velocity dispersion profile and to understand how the 
difference between the observed velocity dispersions and the 
model B predictions arises. If the clouds are decelerated by 
the drag force their orbits become more radially extended 
and the velocity ellipsoids acquire a radial anisotropy. This 
is a general effect and it is true not only in our case, in which 
the clouds initially moved onto circular orbits. So we expect 
(in the region of the gaseous spheroid) the observed velocity 
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Figure 17. Variation of the tangential and radial components of 
the velocity of ionized gas nebulae which are initially supposed to 
move at circular velocity (continuous line) in the potential of the 
Eq. ^ under the effect of the drag force given by the Eq. |s| (a) 
and (b) The dashed lines represent the profiles of the tangential 
and radial components of the nebulae velocity as a function of 
the galactocentric distance after a time t = 2 Myr. (c) The radial 
profile of the velocity tangential component after after t = 0.5, 1 
and 2 Myr. (d) The radial profile of the velocity tangential com- 
ponent after t = 2 Myr for nebulae initially moving at a velocity 
0.75 times the circular one. 



dispersion profile to decrease steeper than the one predicted 
by the isotropic model B. In fact inside 5" Si r 10" the 
F97 ionized gas data show that the gas velocity dispersion 
does not already exceed 50 km s _1 , even if its rotation curve 
falls below the circular velocities inferred from the stellar 
kinematics. Given we do not know the the lifetime and the 
density of the cloudlets we cannot make a definite predic- 
tion. 

The drag effects could explain the differences between 
the observed and predicted velocities, if the decreasing of 
the tangential velocities (showed in Fig. ^) is considered as 
an upper limit. 

A proper luminosity-weighted integration of the tangential 
velocities of nebulae along the line-of-sight has to be taken 
into account for the gaseous spheroid. Being the clouds not 
all settled on a particular plane, the plotted values will be 
reduced by two distinct cosine-like terms depending on the 
two angles fixing the position of any particular cloud. 
Moreover (in the general case) , the clouds are not all starting 
at the local circular velocity and also at any given radius 
we find either 'younger' clouds (just shed from stars) and 
'older' clouds, which are coming from farther regions of the 
gaseous spheroid and are going to be thermalized. Due to 
the nature of the drag force (acting much more efficiently 
on fast-moving particles), it is easy to understand that the 
fast 'younger' clouds will soon leave their harbours, while the 
slow 'older' clouds will spend much more time crossing these 



regions. For instance, in the case of nebulae shed at the local 
circular velocity on a given plane (described by Eqs. and 
pl| and shown Fig [n]) at r = 0.96 kpc (where the circular 
velocity is ~ 230 km s _1 ) we found that quite an half of the 
gas cloudlets are coming from regions between 0.97 and 1.1 
kpc with tangential velocity between 150 km s _1 and 127 
km s _1 . 

This scenario is applicable only outside 3", whether if 
it is applicable also inside the region of the central discrete 
structure revealed by the HST image in Fig. ta, it is not clear. 



5 DISCUSSION AND CONCLUSIONS 

The modeling of the stellar and gas kinematics in NGC 4036 
shows that the observed velocities of the ionized gas, mov- 
ing in the gravitational potential determined from the stel- 
lar kinematics, cannot be explained without taking the gas 
velocity dispersion into account. In the inner regions of 
NGC 4036 the gas is clearly not moving at circular veloc- 
ity. This finding is in agreement with earlier results on other 
disc galaxies (Fillmore et al. 1986; Kent 1988; Kormendy & 
Westpfahl 1989) and ellipticals (CvdM94). 

A much better match to the observed gas kinematics is 
found by assuming the ionized gas as made of collisionless 
clouds in a spheroidal and a disc component for which the 
Jeans Equations can be solved in t he g ravitational potential 
of the stars (i.e., model A in Sect. 



4.1) 



Better agreement with the observed gas kinematics is 
achieved by assuming that the ionized gas emission comes 
from material which has recently been shed from the bulge 
stars (i.e., model B in Sect. |l.l| ). If this gas is heated to the 
virial temperature of the galaxy (ceasing to produce emis- 
sion lines) within a time much shorter than the orbital time, 
it shares the same true kinematics of its parent stars. If this 
is the case, we would observe a different kinematics for ion- 
ized gas and stars due only to their different spatial distri- 
bution. The number of emission line photons produced per 
unit mass of lost gas may depend on the environment and 
therefore varying with the galactocentric distance. Therefore 
the intensity radial profile of the emission lines of the ion- 
ized gas can be different from that of the stellar luminosity. 
The continuum-subtracted Ha+[N n] image of the nucleus 
of NGC 4036 (Fig. ^ confirms that except for the com- 
plex emission structure inside ~ 3" the smoothness of the 
distribution of the emission as we expected for the gaseous 
spheroidal component. 

In conclusion, the 'slowly rising' gas rotation curve in 
the inner region of NGC 4036 can be understood kinemat- 
ically, at least in part. The difference between the circular 
velocity curve (inferred from the stellar kinematics) and the 
rotation curve measured for the ionized gas is substantially 
due to the high velocity dispersion of the gas. 
This kinematical modeling leaves open the questions about 
the physical state (e.g. the lifetime of the emitting clouds) 
and the origin of the dynamically hot gas. We tested the hy- 
pothesis that the ionized gas is located in short-living clouds 
shed by evolved stars (e.g. Mathews 1990) finding a reason- 
able agreement with our observational data. These clouds 
may be ionized by the parent stars, by shocks, or by the 
UV-flux from hot stars (Bertola et al. 1995a). The compar- 
ison with the more recent and detailed data on gas by F97 
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opens wide the possibility for further modeling improvement 
if the drag effects on gaseous cloudlets (due to the diffuse 
interstellar medium) will be taken into account. These argu- 
ments indicate that the dynamically hot gas in NGC 4036 
has an internal origin. This does not exclude the possibility 
for the gaseous disc to be of external origin as discussed for 
SO's by Bertola, Buson & Zeilinger (1992). 

Spectra at higher spectral and spatial resolution are 
needed to understand the structure of the gas inside 3". 
Two-dimensional spectra could further elucidate the nature 
of the gas. 
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APPENDIX A: DATA TABLES 



Table Al. Stellar kinematics along the major axis of NGC 4036 



The stellar (Tab. Al) and ionized gas (Tab. A2) heliocentric 
velocities and velocity dispersions measured along the major 
axis of NGC 4036. 



V 


V 


V 


(j 


OCT 


[ J 


[km s 1 


[km s ] 


firm c — ll 
[Kill h 


[kin s 


\ L ) 


(2) 






{») 


97 Q 
— Z / .y 


1 ftftl 
lOOl 


9*3 
Zo 


1 1 Q 

i iy 


Oft 
Zo 


91 o 


1 ft'30 
lOOU 


1 Q 

iy 


7ft 
1 O 


OO 


17 Q 

— i / .y 


1 ftA1 
lOH 


1 ft 

lo 


on 
yu 


ou 


i ^ ^ 

— 10. o 


i fin*} 


99 
ZZ 


1 ^i 


9 1 
Z 1 


1 1 ^ 

— 1G.O 


1 fi^X^ 
lOOl 


97 
Z i 


ioy 


9ft 
Zo 


19 
— 1Z.U 




9 "1 
Zl 


1 o^ 

lUu 


90 

zy 


1 n q 
— iu.y 


i p.'X'X 

lOOO 


o ■. ) 


1 Q7 

iy / 


9 r i 
Z-O 


Q Q 

— y .y 


1 ^70 
lO I u 


Zo 


lift 

1 lo 


90 

zy 


Q 1 

— y . i 


1 ^09 

ioyz 


97 
Z / 


lOU 


9R 
ZO 


ft A 


1 

lOOl 


99 
ZZ 


"1 ft ft 
1DO 


99 


7 8 


1 ^ft^i 
lOOO 


1 Q 

iy 


1 ftft 
lOo 


9 1 
Z 1 


7 1 


1 ^ftQ 

iooy 


9 3 
ZO 


1 7'\ 
1 I o 


9 ^ 
ZO 


ft A 
— O.l 


1 ^7ft 
10 / o 


OA 
Z 1 


iyo 


9 1 
Z 1 


^ ft 
— O.o 


lOOO 


ID 


1 fift 

lOO 


"1 ("1 

iy 


— O.l 


lO^O 


1 Q 

iy 


"I S^i 
loo 


90 
ZU 


A A 


1 ^ftO 
lOOU 


9ft 
Zo 


90^ 
ZUO 


99 
ZZ 


"X ft 
— o.o 


1 AQft 


"1 7 
1 / 


ZUO 


"1 7 
1 f 


'X 1 
— o.l 


1 ^1 1 
1011 


1 

iy 


1 Q7 

iy / 


"1 

iy 


9 R 
— z.o 


1 ^04 
10U4 


14 


201 


iu 


i s 
— 1 .o 


1 zlftQ 

i^oy 


i 1 

1 1 


1 fti 
lol 


lo 


— 1 1 


1 Aftft 


1 
IU 


"I (\ "I 

1 .) 1 


1 4 


— U. O 


1 AM 
1101 


1 "I 
1 1 


inr: 

iyo 


"1 A 
IT 


("1 


1 Af\7 


i n 

1U 


90Q 
ZU.y 




U.o 


1 ^ftO 


("1 
y 


ZUO 


lo 


l.u 


1 'XFt'x 

lOOO 


11 


lou 


lo 


9 "1 
Z. 1 


1 "XA^ 


1 

iu 


"1 7^ 
1 ( C5 


"1 1 
IT 


9 a 

Z.o 


1 "X'XC\ 
lOOU 


1 9 
1Z 


1 ftQ 

ioy 


lO 


o . ; ) 


1 'XOH 
loZo 


1 'X 
lo 


"I S1 
lol 


ID 


A 1 
'1. 1 


1 ^09 


ID 


"I ("I "I 

i y i 


"1 ft 
lo 


A 8 
T.o 


1 "307 
lou 1 


1 1 

1-1 


1 ft 1 ^ 

1O0 


"1 7 
1 f 


E A 


1 "XClA 


1 1 
11 


in 


90 
ZU 


ft "1 
0. 1 


1 9Q0 

izyu 


9 "1 
Zl 


inc; 
i y ■ ) 


90 
ZU 


ft H 
O.o 


1 9Q0 

izyu 


1 Q 

iy 


1 70 
1 i U 


9 1 
Z 1 


7 A 
i . 1 


1 9ftft 
IZOo 


1 

iy 


1 7^ 
1 I 


90 
ZU 


ft 7 


1 9QQ 

izyy 


1A 
o 1 


907 
ZU 1 


9^ 
ZO 


y . 1 


1 9A1 


O 1 


1 Qft 


9 r i 
ZO 


10 9 
1U.Z 


1 9ftft 
IZOO 


"1 7 
1 i 


110 
1 1U 


9ft 
ZO 


11.2 


1199 


46 


204 


30 


12.4 


1186 


28 


141 


28 


13.7 


1225 


22 


98 


32 


15.2 


1239 


24 


140 


26 


17.0 


1172 


22 


104 


30 


19.0 


1193 


23 


131 


26 


21.1 


1202 


18 


36 


68 


23.7 


1155 


26 


176 


24 


29.1 


1160 


19 


80 


34 



© RAS, MNRAS 000, §-|l7| 



16 Cinzano et al. 



Table A2. Ionized gas kinematics along the major axis of 
NGC 4036 
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